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Recent developments provided evidence that the dimension 2 gluon condensate (A 2 ) is important 
for the nonperturbative regime of Yang- Mills theories (quantized in the Landau gauge). We show 
that it may be relevant for the Dyson- Schwinger approach to QCD. In order that this approach 
leads to a successful hadronic phenomenology, an enhancement of the effective quark-gluon inter- 
action seems to be needed at intermediate (p 2 ~ 0.5 GeV 2 ) momenta. It is shown that the gluon 
condensate (A 2 ) provides such an enhancement. It is also shown that the resulting effective strong 
running coupling leads to the sufficiently strong dynamical chiral symmetry breaking and successful 
phenomenology at least in the light sector of pseudoscalar mesons. 
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1. INTRODUCTION AND SURVEY 

Dyson-Schwingcr (DS) equations provide a prominent 
approach to physics of strong interactions. One of its 
aspects, reviewed Q and exemplified recently by Refs. 
1, 2], consists of ab initio studies of DS equations for 
Green's functions of QCD, typically in the Landau gauge 
(LG) . The other aspect consists of phenomenological DS 
studies (also typically in LG) of hadrons as quark bound 
states, but relying more on modeling; e.g., see a re- 
cent review j3(- Such phenomenological studies have 
mostly been relying on the rainbow-ladder approxima- 
tion (RLA), where generation of dynamical chiral sym- 
metry breaking (D%SB) is well- understood [E S III bl 111 • 
As it has been stressed in, e.g., Refs. IJ0> RLA is the 
leading-order term of a procedure 0, @ that can be sys- 
tematically improved towards less severe truncations of 
DS equations. This general procedure provides a means 
to identify a priori the channels in which RLA is likely to 
work well . Pseudoscalar mesons are the most notable 
among those channels because of the correct chiral QCD 
behavior due to D^SB and Goldstone theorem, since (al- 
most) massless pseudoscalar mesons are reproduced in 
the (vicinity of) chiral limit not only by the exact QCD 
treatment but also by all consistent truncations such as 
RLA [13 . 

Consistent RLA implies Q II II H H that Ansdtze of 
the form 

[K( P t 9 f = -47HW) Df u (p) 7/1 ] e9 [| lu ] hf (1) 

must be used for the interactions between quarks in both 
the gap equation S* _1 = Sq 1 — E for the full quark prop- 
agator S (Sq is the free one) and the Bethe-Salpeter (BS) 
equation for the meson (M) bound-state vertex T M ; i.e., 

PW = J S gh [K] h e f , [r M ] eJ = J [ST*S] gh [K] h J> , (2) 

where, writing schematically, integrations are meant over 
loop momenta while e,f,g,h in Eqs. represent 
spinor, color and flavor indices 0- In LG, the free gluon 
propagator Df v {p) = 8 ab {5 uv - p^p„ /p 2 ) /p 2 . 



In Eq. a c s(p 2 ) is an effective running coupling. 
It is only partially known from the fact that at large 
spacelike momenta (our convention is p 2 > for spacelike 
p), a c jf(p 2 ) must reduce to a pe (p 2 ), the well-known run- 
ning coupling of pcrturbative QCD. However, for p 2 < 1 
GeV 2 , where nonperturbative QCD applies, the inter- 
action is still not known. Thus, in phenomenological 
DS studies, a e ff(p 2 ) must be modeled for p 2 < 1 GeV 2 
- e.g., see Refs. 0, 0, 0, IE 01 ■ There one can see 
that phenomenologically most successful of those mod- 
eled interactions have a rather large bump at intermedi- 
ate momenta, around p 2 ~ 0.5 GeV 2 (e.g., in Fig. Q]see 
a cff(p 2 ) used by Jain and Munczek (JM) |j| and by Maris, 
Roberts and Tandy (MRT) HHQ). In any case, success- 
ful DS phenomenology demands that the modeled part of 
the interaction JIJ be fairly strong; regardless of details 
of the form of the interaction, its integrated strength (for 
p 2 < 1 GeV 2 ) must be fairly high to achieve acceptable 
description of hadrons, notably mass spectra and D%SB 
■ Theoretical explanations of the origin of so strong 
nonperturbative part of the phenomenologically required 
interaction are obviously much needed, cither from the 
ab initio DS studies or from somewhere outside the DS 
approach. This is the main motivation for the present 
paper. 

The ab initio DS studies showed that, in LG, the effects 
of ghosts are crucial for the intermediate-momenta en- 
hancement of the effective quark-gluon interaction. This 
is obvious in the expression for the strong running cou- 
pling a s (p 2 ) in these LG studies 0,05 

a s {p 2 ) = a s {ft 2 )Z(p 2 )G(p 2 ) 2 , (3) 

where a s (fi 2 ) = .9 2 /4tt and Z(fi 2 )G(^ 2 ) 2 = 1 at the 
renormalization point p 2 = /i 2 . The ghost and gluon 
renormalization functions G(p 2 ) and Z(p 2 ) define the full 
ghost propagator D<£(p) = -S ab G{p 2 )/p 2 and the full 
LG gluon propagator Df u {p) = Zip 2 ) D±(pV 

While the ab initio DS studies [J || do find signif- 
icant enhancement of a s (p 2 ), Eq. (jSJ), until recently 
this seemed still not enough to yield a sufficiently strong 
D%SB (e.g., see Sec. 5.3 in Ref. 0) and a successful phe- 
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FIG. 1: The momentum dependence of various strong run- 
ning couplings mentioned in the text. JM Q and MRT 
■3,.6j a e ff(p 2 ) ar e depicted by, respectively, dashed and dash- 
dotted curves. The effective coupling 11811 proposed and ana- 
lyzed in the present paper is depicted by the solid curve, and 
a s (p 2 ) © of Fischer and Alkofer (fit A) by the dotted 
curve. The long-dashed curve is the fit (extrapolated all the 
way to p 2 = 0) of the lattice results of Ref. 



nomenology. Nevertheless, going beyond the ladder trun- 
cation and so getting additional interaction strength from 
dressed vertices, for carefully constructed dressed quark- 
gluon vertex Ansatze, Fischer and Alkofer |2| have re- 
cently obtained good results for constituent quark masses 
and pion decay constant f^, although not simultane- 
ously also for the chiral quark-antiquark (qq) conden- 
sate, which then came out somewhat too large. Thus, 
the overall situation is that there is progress 0, El , but 
that further investigation and elucidation of the origin of 
phenomenologically successful effective interaction ker- 
nels remains one of primary challenges in DS studies Q . 
Here, we will point out that such an interaction kernel for 
DS studies results from combining the DS ideas on «s(p 2 ) 
of the form © Q, and the ideas d Q E El [13 
on the dimension 2 gluon condensate (A^A" 1 ^) = (A 2 ) in 
LG. 



2. CONDENSATES IN GLUON AND GHOST 
PROPAGATORS 

Already a long time ago Refs. [II H HJ El found 
in the operator product expansion (OPE) the (A 2 )- 
contributions to QCD propagators, recently confirmed 
by Kondo ^3- For LG (adopted throughout this pa- 
per), number of colors N c = 3 and space-time dimensions 
D = 4, their results for gluon and ghost propagators 
amount to 



Qg(i/p 2 



2 ° 2 / a2\ 

A = 32 9 ' ' = 



(5) 



(6) 



where tua and ma are, respectively, dynamically gen- 
erated effective gluon (A) and ghost (G) mass. The 
later references also worked out the perturbative 

QCD corrections inducing the logarithmic p 2 -dependcncc 
of these dynamically generated masses, i.e. m A (p 2 ), to 
which we will return and comment on in the next section. 

These now well-established propagator contributions 
©-© then suggest the importance of {A 2 ) for the DS 
approach to hadrons, where propagators, usually in LG, 
are used to get solutions for quark bound states and cal- 
culate observable quantities. Notably, see Ref. ^3 for 
gauge-parameter independent expressions for f n and a 
generalization of the Gell-Mann-Oakes-Renner relation 
(GMOR) that demonstrates gauge-parameter indepen- 
dence of the meson mass. Still, how can (A 2 ) influence 
these observable quantities, since this condensate is not 
gauge-invariant? It turns out 0, E, 0] that in LG, 
(A 2 ) equals a non-local, but gauge-invariant quantity: 
the minimal (with respect to the choice of gauge) value 
of A^A^ 0, integrated over the space-time, indicating that 
(A 2 ) in LG may have a physical meaning. Outside LG, 
besides (A 2 ) other (ghost) condensates of dimension 2 
appear |16| . They very likely cancel the variation which 
(A 2 ) suffers in going to another gauge, since the physics 
behind all these different dimension 2 condensates in dif- 
ferent gauges must be the same: gluon-ghost condensa- 
tion lowers the QCD vacuum energy E, which is a phys- 
ical, gau ge-invariant quantity, to a stable ( " E < 0" ) vac- 
uum 241. 



For g 2 {A 2 ), LG lattice studies of Boucaud et 



Z( P 2 ) = 



Oa(1/p 2 ) 




(4) 



m 

yield the value 2.76 GeV 2 , compatible with the bound re- 
sulting from the discussions of Gubarev et al. 0, El on 
the physical meaning of (^4 2 ) and its importance for con- 
finement. This value gives © tua — 0.845 GeV, which 
will turn out to be a very good initial estimate for niA.G- 
As for the contributions Oa(1/p 2 ) and Og(1/p 2 ) m 
Eqs. © and ©, one expects a prominent role of 
the dimension 4 gluon condensate (F®„F atl ' J ) = (F 2 ), 
which, contrary to (A 2 ), is gauge invariant [25j. Refs. 
[T^. l2(il ] showed that the OPE contributions of dimen- 
sion 4 condensates were far more complicated [2l| than 
found previously |18| : not only many kinds of conden- 
sates contributed to terms oc 1/p 2 , but for many of them 
(gauge-dependent gluon, ghost and mixed ones) there has 
been no assignments of any kind of values yet. Terms 
oc (1/p 2 )" (n > 1) were not considered at all. Thus, at 
this point, the only practical approach is that the con- 
tributions O a (1/p 2 ) and G (l/p 2 ) in Eqs. © and © 
are approximated by the terms oc 1/p 2 and parametrized, 
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where r^(p 2 ), r^ip 2 ), rf(p 2 ), etc., are the terms of mass- 
dimension 0, 2, 4, etc., for the gluon case, and r^(p 2 ), 

Oa(1/p 2 ) ~ — — , Og(1/p 2 ) ~ — — , (7) r 2{v 2 )i r 2{v 2 )i etc., are the terms of mass-dimension 0, 

P 2, 4, etc., for the ghost case. For example, the terms of 

the dimension zero, up to one loop, are 



Ca 
P 2 



where both Ca and Cg are in principle free parameters 
to be fixed by phenomenology. Still, we should mention 
that the effective gluon propagator advocated by Lavelle 
[H would indicate C A ~ (0.640 GeV) 4 for the follow- 
ing reason: for LG and D = 4, the contribution which 
this gluon propagator receives from the so-called "pinch 
diagrams" vanishes, and its [26] Oa(1/p 2 ) contribution 



n 



(F 2 ) 



(P 2 



UN c na s {F 2 } _ (0.640 GeV) 4 
9(N 2 - l)p 2 ~ p 2 



(8) 



stems entirely from the gluon polarization function in 
Ref. |l9j| . provided one invokes some fairly plausible as- 
sumptions, like using equations of motion, to eliminate 
all condensates except (F 2 ). (The quark condensate (qq) 
could also be neglected Hy.) Since Ref. (2?J indicates 
that the true value of a s (F 2 ) is still rather uncertain, 
and since Refs. [2(3, make clear that Lavelle's j2|| 
propagator misses some (unknown) three- and four-gluon 
contributions, we do not attach too much importance to 
the precise value C A = (0.640 GeV) 4 [HEl in Eq. ©, 
but just use it as an inspired initial estimate. Fortu- 
nately, the corresponding variations of Ca still permit 
good phenomenological fits, since we will find below that 
our results are not very sensitive to Ca- 

There is no similar estimate for Cg, but one may sup- 
pose that it would not differ from Ca by orders of magni- 
tude. We thus try C G = C A = (0.640 GeV) 4 as an initial 
guess. It will turn out, a posteriori, that this value of Cq 
leads to a remarkably good fit to phenomenology. 



4(P 2 ) = 



*pc 



(P 2 ) 



' pc 



(^ 2 ) 



«pc(p 2 ) 

a P c(^ 2 ) 



(11) 



(12) 



where 7 and 5 are, respectively, gluon and ghost anoma- 
lous dimensions. The perturbative corrections for the 
Wilson coefficients of dimension 2 have also been calcu- 
lated for the pure Yang-Mills case by Boucaud et al. |2S| 
(for the gluon propagator) and Kondo and collaborators 
[22ll2^ ] (for both the gluon and ghost propagators), and 
since they imply that 



V 



Kp 2 ) 



r£(p 2 



OT 2 ) 



m A(p 2 ) = - m G(p 2 ) 



we write the renormalization functions as 



Z(P 2 ) 
G( P 2 ) 



a P c(p 2 ) 
a P o(M 2 ) 

a P c(p 2 ) 
a pc ([i 2 ) 



™\(P 2 ) , C A (P 2 ) 



+ 



c g (p 2 



(13) 

(14) 
,(15) 



COUPLING ENHANCED BY THE GLUON 
CONDENSATES 



Having set the stage, we are now ready to propose that 
m 2 A = -nig oc (A 2 ) leads to the enhancement of a c g(p 2 ) 
at intermediate p 2 . To derive the running coupling ex- 
hibiting this property, let us first recall the aforemen- 
tioned perturbative corrections to OPE results (0J, © 
and ©. In Eqs. and Q, gluon and ghost renor- 

malization functions are parametrized by the coefficients 
rriA, Ca, and Cg, which are constants at the tree level 
but develop momentum dependence through the pertur- 
bative corrections. To see this, we note that the generic 
forms of the ghost and gluon renormalization functions 
including OPE contributions and perturbative QCD cor- 
rections [22|, 123, |23 can be written as 



Z(P 2 ) 



G( P 2 ) 



1 



r${p 2 ) 



1 



G{„2\ 1 r ?(P 2 ) , r X (P 2 ) 



r^{p 2 ) 



(9) 
(10) 



The perturbative corrections for the Wilson coefficients 
of dimension four and higher have not been calculated 
yet, but we introduced also the notation 



rj(p 2 ) 
ro(P 2 ) 



Ca(p 2 ) 



! (P 2 ) 



r<f(p 2 ) 



= c g (p 2 ) 



(16) 



to point out the correspondence of Eqs. 114|) and i|15|) 
with the relations Q-© and J7J). The latter differ from 
the former just by the absence of the slowly varying log- 
arithmic p 2 -dependence in mA{p 2 ), Ca{p 2 ), Cg(p 2 ) and 
Qpcip 2 ), and of the dots denoting terms of dimension 
larger than four. 

Regarding the perturbatively generated prefactors 
(apc(p 2 )/o!pc(/i 2 ) to the powers of —7 and —S), the forms 
(|14f> and l|15(l are consistent with the corresponding forms 
given by Eqs. (12) and (13) and also (14) in Ref. 
and by Eqs. (41) in Ref. 2]. For N c colors and N f 
quark flavors, the anomalous dimensions of the gluon and 
ghost propagator are respectively given by 7 = ( — 13iV c + 
4N f )/(22N £ - 4N f ) and 5 = -9N C /(UN C - 8N f ) (see, 
e.g., Ref. yQ). This ensures 7 + 25 = —1. The definition 
of the strong running coupling constant, Eq. ||2J), together 
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with Eqs. l(T4^l and thus gives 

1 

a s (p 2 ) = a pe (p 2 ) 



1 + ™ A (p2) + C A (p 2 ) + 



1 p 2 1 p^~ 



(17) 



Neglecting the p 2 -dependence in the coefficients m,A, Ca, 
and Cq, as well as the higher terms in the denominators, 
we finally get 



a s (p 2 ) S3 a pe (p 2 )- 



1 



1 



m A , C A 



(18) 



a eS (p 2 ) = a pc (p 2 ) Z^{p 2 ) G^{p 2 ) 2 



depicted in Fig. \I\by the solid line [for the parameter val- 
ues (J2EJ, discussed below]. The suggestive abbreviations 



and 



Z N P c (p 2 



G Npc (p 2 ) 



C A 



1 



Cg 
P 4 



(19) 



(20) 



for the factors giving the deviation of Eq. I|18|) from 
the perturbative coupling a pc (p 2 ), stress that our ap- 
proximations amount to assuming that nonperturbative 
(Npe) effects are given by the OPE-based results of Refs. 

d 111 E M m 

which in our present case boil down 
to Eqs. Q-©, and by the parametrization J7J. 

Our final expression (|18fl for the running coupling, to 
be used in phenomenological calculations below, does not 
depend on the renormalization scale /z 2 explicitly, al- 
though Z(p 2 ) and G(p 2 ), appearing in the intermediate 
steps, still do. Namely, the same approximation applied 
to the gluon and ghost renormalization functions 114|) 
and l(T5|l . gives 



Z{P 2 ) 
G( P 2 ) 



a P o(M 2 ) 

a P c{p 2 ) 
" P c(^ 2 ) 



1 



C A 



-6 



1 - 



C G 



(21) 
(22) 



plotted in Figs. [21 and as our G(p 2 ) and Z(p 2 ) [again 
for the parameter values (|26|) ]. 

In our figures we also plot some results of lattice and 
ab initio DS studies. Namely, before we turn to our main 
goal, i.e., exploring whether our running coupling l|18l) 
leads to successful phenomenology when used in quark 
gap DS and bound-state BS equations through Eq. 



we will first comment on the comparison of our gluon 
and ghost renormalization functions with some other re- 
sults. In particular, recent Tubingen results of SU(2) 
lattice gauge simulations [eH I34L l35| and of ab initio DS 
studies 0) yQ| agree with each other well for G(p 2 ) and 
quite reasonably for Z(p 2 ). This is seen in Figs. [21 and 
(Also, the results in most recent Tubingen lattice refer- 
ence |35l , which however does not give the corresponding 
fitting formulas, are after proper renormalization some- 
what lower than the results of Ref. ^3 plotted in Figs. 
[21 and [3 and thus agree somewhat better with the dis- 
played ab initio DS results 0, H3-) As can be seen in 
Fig. [3J there are some other lattice results 0, 13^ |3^| 
for the gluon renormalization function which agree with 
Z(p 2 ) from DS approach 0, H(j even better, but they 
do not give the ghost renormalization function. Thus, 
Tubingen lattice results 0, 0, [3f| are presently of par- 
ticular interest, because they give both Z(p 2 ) and G(p 2 ). 
Admittedly, there is a caveat: while the ab initio DS 
studies 0, E3 do not have problems with reaching low 
momenta and in fact make strong statements about the 
asymptotic behavior in the p 2 — > limit, the lattice data 
1 121] do not reach very low momenta. The lowest data 
point for G{p) as well as a{p) in Ref. ^3 is at p 2 ~ 0.36 
GeV 2 , so that one must keep in mind that for lower p 2 , 
the corresponding lattice-data-fitting curves in Figs. [21 
and [3J and therefore also in Fig. are just extrapo- 
lations. Nevertheless, presently most important is that 
comparing the long-dashed and dotted curves in Fig. ^ 
shows that the respective running couplings © follow- 
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FIG. 2: The momentum dependence of the ghost renormal- 
ization function G(p 2 ). The solid curve is our result 121'} . The 
dotted curve is a result of the ab initio DS study 01 ■ Con- 
cretely, it depicts the fit A in Eq. (41) of Ref. p| evaluated 
at u — 5 GeV, for which value the dotted curve agrees rather 
well (in the displayed momentum range) with the long-dashed 
curve, representing Eq. (4) of the lattice Ref. [T^l . renormal- 
ized at u = 5 GeV. Our result 12211 is therefore plotted for the 
same value of u. 
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FIG. 3: The momentum dependence of the gluon renormal- 
ization function Z(p 2 ). The solid curve is our result 1 211 . 
The densely-dotted curve is the ab initio DS result (fit A) 
from Refs. |H30|. As in Fig. |U both the ab initio DS result 
and our Eq. 12 H are plotted for fi = 5 GeV. The other curves 
pertain to some recent lattice results [T5 . l3ll 135, 133^ . The 
ones that agree with the a6 initio DS result |2,|3(| very well, 
are those of Leinweber et al. [3l| (their Eq. (5.14), depicted 
here by the sparsely-dotted curve), and Eq. (3) of Iida et al. 
[3^1 (dash-dotted curve) fitting the lattice data of Refs. [33| . 
The third lattice result for Z(p 2 ) is displayed by the long- 
dashed curve. It corresponds to Eq. (3) of Ref. 12], which 
is the fit (extrapolated all the way to p 2 = 0) to the lattice 
data of Ref. |T^|. For all lattice results, the renormalization 
condition Z(pL 2 ) — 1 is imposed at /i = 5 GeV. 



ing fr om these "lattice" 0, 134I l35| and 11 ab initio DS" 
[2l l30j renormalization functions typically do not differ by 
more than the factor of two. Thus, in spite of the men- 
tioned caveat, these lattice results 0, 0, 111] and the 
aforementioned D^SB scenario of Fischer and Alkofer 2] 
support each other. Still, the behavior of QCD propaga- 
tors, especially the ghost ones, and the resulting running 
coupling, is not a closed issue yet, so that the presently 
proposed scenario should also be considered although it 
is not supported by lattice results. 

The examples |M 111 HJ^lof lattice results differing 
from the Tubingen lattice 12, 35] and ab initio DS results 
@,|3(j are not only the relatively old ones such as those of 
Suman and Schilling on the ghost renormalization func- 
tion (which abruptly falls for the very smallest probed 
momenta, possibly indicating the infrared vanishing be- 
havior) |3a, l3j. but also some of the most recent ones, 
such as Ref. |2g, where Landau gauge lattice calcula- 
tions give the strong running coupling which, supposedly 
due to instanton effects, decreases 1 at small momentum 



roughly as p |38|, |3_ 

Some other, quite independent methods, also give the 
QCD running coupling vanishing at small p 2 |40j al- 
though not so fast as our form (|18|) . Now, we want to 
make clear that we do not argue that our results are 
another indication that the running coupling indeed van- 
ishes as p 2 — > 0, because we are aware (as we comment 
in more detail below) that the behavior of our running 
coupling at very small p 2 is just an artifact of the way 
Eq. (I18|) was derived. Fortunately however, it turns 
out that the small-p 2 behavior does not influence much 
our final, observable results (in contradistinction to in- 
termediate p 2 's). To see all this, let us discuss in de- 
tail the behavior of our form i|18|) and especially the 
possible objections to it. The first, less serious one is 
that a pe {p 2 ) ultimately hits the Landau pole as p 2 gets 
lower. However, we handle this as in other phenomcno- 
logical DS studies 0, 0, IE E3, E3, E3, E3, E^] , where this 
pole is shifted to timelike momenta in all logarithms: 
ln(p 2 / A QCD) M x o + P 2 /Aq CD ). (Dynamical gluon 
mass can provide the physical reason for this |46f : i.e., 



xq oc m A 



Mqcd 



10.) For a pe (p ) we use the two- 
loop expression used before by JM H and our earlier DS 
studies, e.g., Refs. [H |H El EH This means the 
infrared (IR) regulator (to which all results are almost 
totally insensitive) is xo = 10, and Aqcd = 0.228 GeV. 
The parameters of a pe (p 2 ) are thereby fixed and do not 
belong among variable parameters such as C a and Cq. 

Back to the possible objections: the second, in 
the present context the more serious one is that we 
cannot in advance give an argument that the factor 
Z N P°(p 2 )G Npo (p 2 ) 2 in the proposed a cS (p 2 ) © indeed 
approximates well nonperturbative contributions at low 
p 2 (say, p 2 < 1 GeV 2 ), but can only hope that our results 
to be calculated will provide an a posteriori justification 
for using it as low as p 2 ~ 0.3 GeV 2 [since Eq. (|T%|) 
takes appreciable values down to about p 2 ~ 0.3 GeV 2 ]. 
Of course, ^ Npe (p 2 ) and G Npe (p 2 ) must be wrong in the 
limit y 2 — > 0, as the OPE-based results (QJ-lO of Refs. 
El El 13 IM E3 certainly fail in that limit. Thus, 
the extreme suppression for small p 2 is an unrealistic 
artifact of the proposed form H18|) when applied down 
to the p 2 — * limit. Nevertheless, because of the inte- 
gration measure in the integral equations in DS calcu- 
lations, integrands at these small p 2 do not contribute 
much, at least not to the quantities (such as (qq) con- 
densate, meson masses, decay constants and amplitudes) 
calculated in phcnomcnological DS analyses. Hence, the 
form of a e ff(p 2 ) at p 2 close to zero is not very impor- 
tant for the outcome of these phenomenological DS cal- 
culations. 2 This is because the most important for the 



1 Their strong running coupling becomes roughly 0.1 or smaller at 
p 2 ~ 0.16 GeV 2 , below which p 2 the lattice evaluation was found 



unreliable l38l . 

2 This is supported by, e.g., Fischer l47f . who found that nearly 
all dynamically generated mass is produced by the integration 
strength above p 2 = 0.25 GeV 2 . 
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success of phenomenological DS calculations seems the 
enhancement at somewhat higher values of p 2 - e.g., see 
the humps atp 2 ~ 0.4 to 0.6 GeV 2 in the JM Q or MKT 
interaction in Fig. Our a e ff(p 2 ) lfT%l) . the solid 

curve in Fig. ^ exhibits such an enhancement centered 



around p 



t /2. This enhancement is readily under- 



stood when one notices that Eq. I|18|) has four poles, 

(P 2 )i,2 - \ ( mi ± i y/4C G -m\ ) , (23) 



(P 2 ) 3 ,4 = - ( -mi ± i J4C A - m 4 ) , (24) 



in the complex p 2 plane. For min{Cc Ca} > m\/A there 
are no poles on the real axis, but saddles between two 
complex conjugated poles. For the DS studies, which are 
almost exclusively carried out in Euclidean space, space- 
like p 2 , i.e., p 2 > is the relevant domain and is thus 
pictured in Fig. There, the maximum of a c s(p 2 ) i|18|) 
at the real axis is at p 2 w mi/2, i.e., the real part of its 
double poles (p 2 )i,2 coming from G Npe (p 2 ) 2 . The height 
and the width of the peak is influenced by both Cq and 
m,A- The enhancement of a e s(p 2 ) (|18|l is thus determined 
by (A 2 ) through Eq. ©, and by the manner this conden- 
sate contributes to the ghost renormalization function. 

We are thus motivated to use this form (|18|) of a c g(p 2 ) 
for all p 2 to test its success in DS calculations. We are 
aware of the shortcomings due to its oversimplified char- 
acter, but its study helps to answer whether the (A 2 ) 
condensate, which has recently attracted so much atten- 
tion, may be important for the enhancement of the effec- 
tive interaction needed for successful DS phenomenology. 
The results presented below indicate that the {A 2 ) con- 
densate may indeed provide an important mechanism not 
considered so far. 




2.5 
p 2 [GeV 2 ] 



FIG. 4: The momentum dependence of the dynamically gen- 
erated quark mass M(p 2 ) for u and d quarks. The solid curve 
is our result for the parameters giving the second line of Ta- 
ble I, but our M(p 2 ) depends in fact very little on the small 
explicit chiral symmetry breaking mass parameters fh u and 
fhd of the very light u and d quarks as long as their values 
are at all realistic. The dotted curve is the ab initio DS result 
HIH. The short-dashed curve is the M(p 2 )-fit of Ref. @ 
to the extrapolation of their lattice data to the chiral limit. 
The dash-dotted curve is the similar result from another lat- 
tice calculation, namely the fit of M(p 2 ) from Fig. 14 of Ref. 

El. 



(a few %) modifications to provide a very good descrip- 
tion of the light pseudoscalar sector. Concretely, now we 
will [both in the chiral limit (\ um ) and realistically away 
from it] quote only results obtained for the parameter set 



4. END RESULTS WITH DISCUSSION AND 
CONCLUSION 

We solved the gap and BS equations J5J) for quark pro- 
pagators 

m i 7 -p„4(p 2 ) + £(p 2 ) ~ i 7 • p + Af(p 2 ) 1 ' 

and for pseudoscalar meson qq (q = u, d, s) bound-state 
vertices (r M ) in the same way as in our previous phe- 
nomenological DS studies, e.g., Refs. [il l 13. 14H H l45j . 
This essentially means as in the JM approach [j], ex- 
cept that instead of JM's a e ff(p 2 ), Eq. (|18|l is em- 
ployed in the RLA interaction £[]). We can thus im- 
mediately present the results because we can refer to 
Refs. 0, EJ IH E! EH for all calculational details, 
such as procedures for solving DS and BS equations, all 
model details, as well as explicit expressions for calcu- 
lated quantities and inputs such as the aforementioned 
IR-regularized a pe (p ). 

These calculations show that the initial rriA , C'a and 
Cg estimates motivated and given above, need only slight 



C A = (0.6060 GeV) 4 = C G , m A = 0.8402 GeV , (26) 

while a broader investigation of parameter dependence 
shows the following, i) The results are only weakly sensi- 
tive to moderate variations (up to the factors of 2 to 1 /2) 
of Ca- H) Contrary to that, the results are very sensi- 
tive to mA and Cg, since they determine the peak of 
our a cS (p 2 ) OH). However, between Cg in - (0.6 GeV) 4 
and C£? ax ~ (0.9 GeV) 4 there are many pairs of these 
quantities which give fits comparable (within a percent) 
to that resulting from the values l|26|) . as long as they 
approximately satisfy the linear relation 

(C G ) 1/4 = 0.7742 m A ~ 0.0442 GeV . (27) 

Thus, the two parameters ruling the strength of a c s{p 2 ) 
are not independent. 

Already the chiral-limit results are very satisfactory: 
f 7t ± = f n o = f„ = 90.5 MeV, the ir° — ► 77 chiral-limit 
amplitude T^( X Hm) = 1/(4^/^) = 0.280 GeV" 1 and 
the chiral condensate (qq) = (—217 MeV) 3 . We also get 
the correct QCD chiral-limit behavior: massless qq pseu- 
doscalars and satisfied (within ~ 4%) GMOR. 
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In the chiral limit, where the quark mass is purely dy- 
namically generated since the bare (and current) quark 
masses vanish, the only parameters are tua, Ca and Cq- 
Away from the chiral limit, chiral symmetry is explicitly 
broken by the nonvanishing bare mass parameters rh q of 
light quarks (q = u, d, s) entering the quark-propagator 
gap equation and the qq BS equation. For the very light 
quarks u and d, the dynamically generated quark masses 
of u- and e?-quarks away from the chiral limit are practi- 
cally the same as in the chiral limit. In Fig. 0] the solid 
curve depicts our results for the momentum dependence 
of these dynamical masses M(p 2 ) = B(p 2 ) / A(p 2 ) of u- 
and d-quarks (in the isosymmetric limit). Fig. ^ also 
presents the lattice results for the u- and (i-quark dy- 
namical masses of Refs. which give their fits to 
the lattice data points in terms of analytic, closed-form 
expressions. Our results for the dynamical quark mass 
M(p 2 ) agree very well with both the ab initio DS results 
[lEfflf and the lattice results jSE^. 

As the first attempt to depart realistically from the 
chiral limit, we adopt without change the fh q values ob- 
tained earlier by JM 4] in a very broad DS fit of the me- 
son phenomenology (with their a e g |4|), i.e., fh u = rhd = 
3.1 MeV and rh s = 73 MeV. Already the corresponding 
results for the masses and decay constants of pions and 
kaons and the 7r° — > 77 decay amplitude T^q ', given in the 
first line of Tabled show a very good agreement with ex- 
perimental values, except for the kaon mass. The second 
line shows that just a slight re-adjustment of the quark 
masses, to fh u = rhd = 3.046 MeV, fh s — 67.70 MeV, is 
enough to get an almost perfect fit to the pion and kaon 
masses. 



The presented results allow us to conclude that the di- 
mension 2 gluon condensate (A 2 ) provides an enhanced 
effective interaction a c g(p 2 ) which leads to a sufficiently 
strong D^SB, pions and kaons as (quasi-)Goldstonc 
bosons of QCD, and successful DS phenomenology at 
least in the light sector of pseudoscalar mesons. This 
opens the possibility that instead of modeling a c s(jp 2 ), 
its enhancement at intermediate p 2 may be understood 
in terms of gluon condensates. 







M K+ 


fK+ 


137 [GeV- 1 ] 


136.17 


93.0 


516.28 


112.5 


0.256 


134.96 


92.9 


494.92 


111.5 


0.256 


134.98 


92.4 ±0.3 


493.68 


113.0 ± 1.0 


0.274 ±0.010 



TABLE I: The masses and decay constants of pions and 
kaons and ir° — > 77 decay amplitude Tjo\ for the param- 
eter values 1261 and our a e g(p 2 ) 1181 . The JM Q quark 
bare masses fh u = rhd = 3.1 MeV, m s = 73 MeV give 
the first line. The second line are the results obtained with 
rhu = rhd = 3.046 MeV , fh s = 67.70 MeV. The last line gives 
the experimental values. (The distinction between neutral 
and charged mesons applies only to this line, as we calculate 
in the isosymmetric limit.) Everything is in MeV except T^o ■ 
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